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An investigation on the effect of oil on the onset of nucleate pool boiling of R-124 
from a single horizontal tube was conducted at saturation temperature of 2.2 °C. Pure R- 
124 and R-124/oil mixtures of 3% and 10% (by weight) miscible alkylbenzene oil were 
used. The tubes tested were: (1) smooth tube, (2) 19 fins per inch and (3) 26 fins per inch 


low-integral finned (GEWA-K) tubes, and (4) porous-coated (HIGH FLUX) tube. The 





effect of dissolved gases, subcooling and pressure on the onset of nucleate boiling were 
also investigated. 

An oil concentration of 3% tends to delay the onset of nucleate pool boiling 
(compared to pure R-124) on the smooth tube and GEWA-K 19 fins per inch tube, but not 
on the GEWA-K 26 fins per inch and HIGH FLUX tubes. The reason for this is not known 
precisely but is very repeatable. A 10% oil concentration tends to delay the onset of 
nucleate pool boiling on all tubes tested. This is due mainly to the increase of surface 
tension, and saturation temperature with increase in oil concentration. The presence of 
dissolved gases in the pool tends to lower the onset of nucleate pool boiling by increasing 


the number of entrapped vapor nuclei. Pressurization and subcooling tend to increase the 





onset of nucleate pool boiling by deactivating potential nucleation sites. 
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I. INTRODUCTION 


A. BACKGROUND 

CFCs are manmade chemicals of chlorine, fluorine and 
carbon. They are used today in numerous applications, such as 
refrigerants for air-conditioning systems, aerosol 
propellants, fire retardants, foam blowing agents, degreasing 
solvents and for defluxing printed circuit boards [Ref. 1}. 

These fluids are generally inert, very stable and do not 
react with other elements in the lower atmosphere. However, 
their stability allows them to reach the stratosphere where 


the earth’s protective ozone layer is located. Here the CFC 


molecules are broken apart by intense ultraviolet radiation, 


producing free chlorine atoms that react with ozone molecules, 
converting them to ordinary oxygen atoms. Since a chlorine 
atom survives the ozone-destroying reaction, it is therefore 
free to repeat the process which can last for a long time - 
one CFC molecule being capable of destroying up to 100,000 
ozone molecules [Ref. 1]. 

Recent studies indicate a continuous reduction of the 
ozone layer worldwide ovez the last 17 years causing an 
increase of ultraviolet radiation with very short wavelengths 
(called UV-B). A long-term increase in UV-B radiation will 


increase the risks to human health, reduce plant growth, cause 








damage to the marine environment and increase atmospheric 
abundance of chemically reactive compounds, thus aggravating 
environmental pollution problems [Ref. 2]. 

With the global concern of the depletion of the ozone 
layer, an international agreement (called the “Montreal 
Protocol" [{Ref. 2]) was reached in Sept 1987. The Protocol 
calls for the freeze in production and consumption of ozone- 
depleting CFCs (to 1986 levels starting in July 1989; this is 
followed by further cutbacks in production and consumption to 
a 20% reduction by 1993, anda further 50% reduction in 1998). 
All CFCs will be completely phased out in 2000 as agreed upon 
by UNEP (United Nations Environment Program) members ina 1990 
progress meeting held in London [Ref. 3]. 

These phase-out dates were revised in November 1992, when 
the Protocol signatories met again in Copenhagen. This latest 
meeting called for 75% reduction of CFCs by 1994 and complete 
phaseout by 1996 (except for essential uses) [Ref. 4]. 

The phaseout of CFCs and particularly CFC-114 (R-114) is 
a major concern of the United States Navy because of the 
commitment the Navy has made to its use on existing ships and 
submarines [Ref. 5]. In compliance with the Montreal 
Protocol, and U. S. Legislation, the Naval Surface Warfare 
Center (NSWC) is currently pursuing research on alternative 
refrigerants. There are two principal alternatives currently 
developed that the Navy is considering to replace existing 


refrigerants: (1) hydrofluorocarbons (HFCs) like R-134a is the 








planned replacement for R-12, since R-134a has zero ozone 
depletion potential (ODP), and (2) hydrochlorofluorocarbons 
(HCFCS) and in particular, HCFC-124 (R-124) which has been 
identified as an intermediate "drop-in" replacement for R-114. 

The Naval Postgraduate School in support of NSWC started 
investigating the nucleate pool boiling characteristics of 
HCFC-124 (R-124) in 1992, Bertsch [Ref. 6] evaluated and 
compared the heat transfer performance of R-124 to R-114 using 
various refrigerant/oil mixtures on a smooth tube and four 
enhanced tubes. He reported that: (a) significant heat 
transfer improvement (>50%) was obtained in pure R-124 
compared to pure R-114, (b) the largest enhancement over all 
heat flux ranges was obtained with the HIGH FLUX tube for pure 
R-124, (c) enhancement in heat transfer performance was 
realized from the smooth and finned tubes with increasing oil 
concentration in R-124 up to 6%. He also investigated the 
effect of oil on the onset of nucleate boiling from a smooth 
tube and a HIGH FLUX tube and found that addition of oil 
delays the point of incipience for the smooth tube, but no 
measurable effect was observed for the HIGH FLUX tube. 

The onset of nucleate pool boiling is a very unstable 
process as pointed out by You et al. [Ref. 7], Memory and 
Marto [Ref. 8], and Bertsch [Ref. 6]. The purpose of this 
thesis is to continue the work with R-124 and to study the 


onset of nucleate boiling under variovs conditions. 


B. OBJECTIVES 

The objectives of this thesis are as follows: 

1. Investigate in greater detail the effect of oil on the 
onset of nucleate pool boiling of R-124 from a smooth and 
enhanced single horizontal tube. 

2. Investigate the effect of dissolved gases, subcooling 
and pressure on the onset of nucleate pool boiling of R-124 


from a single smooth and enhanced horizontal tube. 











II. LITERATURE SURVEY 


A. POOL BOILING CURVE 

Pool boiling occurs when boiling takes place from a heated 
surface immersed in a large pool of quiescent liquid. It is 
used in a variety of applications such as metallurgical 
quenching, flooded shell-and-tube evaporators, and immersion 
cooling of electronic components just to name a few. Figure 
2.1 from Bar-Cohen [Ref. 9] shows the well-known pool boiling 
curve, and defines the various heat transfer regimes. 

1. Single Phase Natural Convection 

This regime (a-b) is in effect as wall superheat is 

increased from an initial low value until just before 
nucleation occurs and the first bubbles form on the heated 
surface. It should be noted that no bubbles are formed during 
natural convection, and that heat is transferred due to the 
buoyancy-driven convection of the liquid. 

2. Boiling Incipience 

Once nucleate boiling starts (point b) and bubbles 

begin to grow and depart from various nucleation sites on the 
heated surface (as depicted in Figure 2.2 from Thome [Ref. 
10]), onset of nucleate boiling (ONB) or incipient boiling has 
been reached. According to Bar-Cohen [Ref. 9], shortly after 


incipience, nucleation frequently propagates rapidly across 
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the heated surface and reduces the wall superheat to the vaiue 
associated with established nucleate boiling (point c). This 
rapid decrease in wall superheat has been described as thermal 
overshoot, or superheat excursion. 

Incipient superheat excursion was first reported by Corty 
and Foust [Ref. 11], and Bankoff et al. [Ref. 12], ain the 
boiling of organic fluids. Bergles and Chyu [Ref. 13] 
compared the nucleate pool boiling characteristics of a HIGH 
FLUX tube and a smooth tube in distilled water and R-113. 
Their results showed that both the HIGH FLUX and smooth tube 
displayed a significant temperature overshoot at ONB. Marto 
and Lepere [Ref. 14] in their investigation of the pool 
boiling of R-113 and FC-72 from smooth and enhanced tubes 
reported that the temperature overshoot was less severe in FC- 
72 than in R-113. Memory and Marto [Ref. 8] investigated the 
effect of-oil on the incipience of nucleate boiling for a 
smooth tube and six enhanced tubes (low finned and modified 
integral-fin) and found that in the natural convection region, 
all seven tubes exhibited similar behavior with and without 
oil. They also found that oil concentration seemed to have no 
systematic effect on ONB for the smooth and low-integral fin 
tubes, and that the ONB occurred at lower heat fluxes for the 
structured surfaces. Bertsch [Ref. 6] investigated the effect 
of oil on the onset of nucleate boiling of R-124 from a smooth 
tube and a HIGH FLUX tube and found that addition of oil 


delays the point of incipience for the smooth tube, but found 








no measurable effect for the HIGH FLUX tube. Several other 


researchers have investigated ONB (see Bar-Cohen and Simon 
[Ref.15]). Onset of nucleate boiling is a very random process 
(You et al. [Ref. 7])). It is a:tected by many factors 
including surface roughness, pressure-temperature history, 
contact angle, dissolved gases, subcooling, fluid properties 
(summarized by Brauer and Mayinger [Ref. 16]), and oil 
addition. All these factors will be discussed in more detail 
later in this chapter. 
3. Nucleate Boiling Region 

As heat flux is further increased from point c, larger 
patches and slugs of vapor depart from the heated surface (see 
Figure 2.2). This region (curve c-d) is commonly called the 
fully developed nucleate pool boiling region. Heat can leave 
the surface as sensible heat in a superheated liquid or as 
latent heat in vapor bubbles. Sustained improvement in heat 
transfer is realized throughout this region. As heat flux is 
further increased, a maximum value of heat flux called the 
Critical Heat Flux (CHF) or Departure from Nucleate Boiling 
(DNB) is attained (point d). 

4. Transition Boiling 

Increase in wall superheat beyond the critical heat 
flux results in a decrease in heat flux. The resulting region 
(curve d-g) is called the transition boiling region. 


According to Thome (Ref. 10], this regime is characterized by 





partial and intermittent contact of the liquid with the heated 


surface, which sign ficantly reduces the heat transfer from 
the heated surface to the fluid. Point g 1s calied the 
Minimum heat flux, and at this wall superheat, the liquid is 
no longer in contact with the heated surface because of the 
rapid rate of vapor generation in the region next to) the 
surface. 
5. Film Boiling 

As wall superheat is increased further from the 
Minimum heat flux (point g), the film boiling region is 
attained. This region is characterized by a stable vapor film 


according to Thome (Ref. 10] that blankets the heated surface 


and separates the liquid from the wall. This region is 
OD O Vapor 
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Figure 2.2 Nucleate Pool Boiling Regimes 


normally avoided since the melting temperature of the heat 


source can be attained quickly. 


B. BOILING NUCLEATION THEORY 

Corty and Foust (Ref. 11], and Bankoff [Ref. 12] were the 
first to relate the presence of vapor bubbles with 
microcavities on the heated surface. They postulated that 
boiling is initiated from these pre-existing trapped vapor 
embryos. Bankoff [Ret. 12] showed that the free energy, or 
work of formation required to create a vapor nucleus at a 
planar solid surface in a liquid (heterogeneous nucleation) 
was less than that required for homogeneous nucleation in the 


liquid phase alone. 


C. VAPOR TRAPPING MECHANISM 
The process of trapping vapor or gas was first considered 


by Bankoff (Ref. 12]. He hypothesized that during the flow of 


liquid over a surface with scratches and cavities, air and/or 


vapor could be trapped in the wedge-shaped grooves if the 
contact angle of the flowing liquid on the solid surface was 
greater than the wedge angle as illustrated in Figure 2.3 
provided by Carey [Ref. 17]. Referring to Figure 2.3(a), 
Carey [Ref. 17] stated: 
When a liquid film flows over a gas-filled cavity, the 
contact angle with the downslope tends to be 
maintained as the liquid begins to fill the cavity. 


Since the liquid front is convex, the "nose" of the 
liquid front will reach the opposite wall before the 





contact line reaches the bottom of the cavity if the 
contact angle is greater than the groove angle 2y. 


The condition for entrapment of gas by the advancing 


liquid front can therefore be stated as: 
6,>2y (2.1) 


where @, is the advancing contact angle that may be 
Significantly greater than the static or equilibrium contact 
angle. 

Now, considering Figure 2.3(b) for a receding liquid 


front, liquid entrapment in the cavity may be stated as 


6,<180°-2y (2.2) 


where @® is the receding contact angle. The difference 


between the advancing contact angle and receding contact 


+— 
Liquid 


Q,, «180°-27 
90“-Y 


Contact Line 
Groove Angle * 2Y 


(a) (b) 


Figure 2.3 Vapor and Liquid Entrapment Schematic in a 
Cavity of (a) an advancing and (b) receding liquid film 








angle during repeated wetting of the cavity due to an 
advancing and receding liquid film is called “contact angle 
hysteresis". Carey [Ref. 17] stated that this hysteresis 
effect may result in 0, being much lower than the static or 
advancing contact angles, and may therefore serve to increase 
the tendency to trap vapor in surface cavities. 

Bankoff [Ref. 12] noted that cavities can be classified 
into one of four categories: (1) cavities that trap vapor/gas 
only, (2) cavities that trap liquid only, (3) cavities that 
trap both liquid and vapor/gas, and (4) cavities that do not 
trap liquid, vapor or gas. Carey [Ref. 17] said that category 
of a cavity can be based on whether equations 2.1 and 2.2 are 
satisfied, and that contact-angle hysteresis enhances the 
probability that cavities that trap gas (category 1 or 3) will 


be present on the solid surface. 


D. NUCLEATE BOILING HEAT TRANSFER MECHANISM 

According to Chongrungreong and Sauer [Ref. 18], nucleate 
boiling starts when the temperature of the heated surface 
exceeds the saturation temperature of the liquid by a few 
degrees. Adjacent to the surface, a thin layer of superheated 
liquid is formed in which bubbles nucleate and grow. 

1. Boiling On Plain Surfaces 

Hsu and Graham [Ref. 19] investigated extensively the 

physical processes involved in nucleate pool boiling from 


plain surfaces and identified the following principal heat 
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transfer mechanisms as depicted in Figure 2.4 from Thome 
(Ref. 20]: 
a. Bubble Agitation 
As displayed in Figure 2.4, the motion imparted to 
the liquid from the growth and departure of bubbles improves 
the liquid-phase convection. Thome [Ref. 20] stated that a 
forced convection process is created in which the velocity of 
the liquid near the heated wall is the same or close to the 
velocity of the growing bubble, and that heat is transported 
as sensible heat in the superheated liquid. 
b. Thermal Boundary-Layer Stripping 
Thome [Ref. 20) said that hydrodynamic drag of the 
departing bubbles periodically scrubs the thermal boundary 
layer formed by transient heat conduction into the surrounding 


liquid. The sensible heat transport rate is controlled by the 


quantity of superheated liquid removed from the thermal 


boundary layer. 
c. Evaporation 

A vapor bubble growing on the heated wall at a 
nucleation site is formed by the vaporization of superheated 
liquid surrounding the bubble and by thin film evaporation in 
the microlayer trapped beneath the bubble. Heat is removed 
away as latent heat. 

The actual boiling process is a combination of 


these three heat transfer mechanisms. According to Thome 
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Pigure 2.4 Nucleate Boiling Mechanisms 


(Ref. 20], the heat removed, however, is not a simple 
summation of these three individual contributions since they 
interact and “overlap” in the heat carried away. 
2. Boiling on Enhanced Surfaces 

Surface geometry can affect significantly the heat 
transfer mechar:isms. As before, heat can leave the surface as 
latent heat i: vapor/bubbles or as sensible heat in 
superheated liquid. With enhanced surfaces, Thome [Ref. 21] 


said that latent heat can be transferred both inside reentrant 
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channels and by further growth of the vapor bubbles on the 
outside surface. Sensible heat can be removed by the liquid 
flowing into and out of the reentrant channels as well as on 
the external surfaces of the enhancement. 
a. Evaporation 

According to Thome [Ref. 20], this process is 
composed of three parts as depicted in Figure 2.5: (1) thin 
film evaporation resulting from the conduction and/or 
convection of heat across the thin liquid film, (2) capillary 
evaporation which results from the vaporization of liquid at 
menisci in the reentrant cavities, and (3) external 
evaporation in which the bubble continues to grow in the 
thermal boundary layer at the outer edge of the enhancement. 

b. Convection 

This process is influenced by: (1) the two-phase 
flow boiling process (first proposed by Stephan and Mitrovic 
{[Ref. 22]) inside the reentrant channels where bubbles grow, 
circulate and eventually squeeze out to the surface and liquid 
fills in, (2) the external convection process which is 
controlled by bubble agitation and the vapor-liquid exchange 
mechanism as with plain surfaces (D.1) of Thome [Ref. 21]. 
Much analysis remains to be done on the mechanisms affecting 


overall heat transfer when boiling from enhanced surfaces. 
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Figure 2.5 Evaporation Process on an Enhanced 
Surface 


EZ. FACTORS AFFECTING INCIPIZENT BOILING SUPERHEAT 
1. Dissolved gases 
Bar-Cohen [Ref. 9] noted that the formation of a 
vapor/gas embryo in a surface cavity creates a void into which 
superheated liquid can evaporate and into which dissolved 
gases can flow. Volume of the gaseous mixture increases 
until pressure equilibrium with the surrounding liquid is 
attained by each bubble. 
Assuming the bubble to be hemispherical and 
recognizing the possible presence of dissolved gases, a force 


balance at the bubble interface for steady or quasi-steady- 


state conditions can be stated ([{Roshenow [Ref. 23]) as: 
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p,- P= 228 -p, (2.3) 


Using the Clausius-Clapeyron relation, Roshenow (Ref. 23] 
further derived an approximate relation for the incipient 


superheat. for small pressure differences, (20/r << p) as: 


TsatVrgi/ 20 
= = sat” f - (2.4) 
(Ty Teaet i Deg ( r P-P» P,) 


Murphy and Bergles [Ref. 24] observed that the 
presence of gas in an embryonic bubble trapped in a surface 
cavity or floating in the liquid can substantially reduce the 
vapor pressure needed to preserve the integrity of the bubble 
thus lowering the incipient superheat. They found that for R- 
113, which was well-mixed with air at room temperature and 1 
bar absolute pressure, the saturation temperature was only 27 
°C versus the normal 48 °C for the pure substance. 

According to You et al. [Ref. 7], however, the 
presence of dissolved gases does not necessarily guarantee a 
reduction in the incipient superheat for nucleate pool 
boiling. In their experiments using FC-72 boiling from a 
cylindrical, sputtered-silicon heater, the incipient superheat 
for degassed liquid varied from 18 °C to 38 °C. 

2. Surface Wettability 

The incipience mechanism of a highly wetting liquid 

was investigated by Bar-Cohen and Simon [Ref. 15], Marsh and 


Mudawar [Ref. 25], and Tong et al. [Ref. 26]. They found the 








following significant differences between highly-wetting 
liquids (such as refrigerants and dielectrics) and a poorly- 
wetting liquid (such as water). With highly-wetting liquids 
there will be fewer active embryos (gas retained in cavities) 
on a given surface than with poorly-wetting liquids. This is 
due to the increased ability of fluids with small wetting 
angle to effectively flood the cavities. Another difference 
is that the residual bubble in highly wetting liquids is small 
and therefore requires a larger superheat for its initial 
growth. 

According to Tong et al. [Ref. 26], parameters such as 
volume of trapped vapor/air in a cavity, critical bubble 
radius, incipient superheat, and the superheat excursion at 
the onset of nucleate boiling are strongly affected by surface 
wettability. Brauer and Mayinger [Ref. 16] state that 
increase in surface wettability results in a decrease in the 
vapor trapping capability of the cavities on a surface, which 
in turn delays the onset of nucleate boiling. 

Highly-wetting liquids like chorofluorocarbons have 
been observed by several investigators to have high wall 
superheats at the onset of nucleate boiling. Marto and Lepere 
[Ref. 14] reported that due to the high wettability of R-113 
on copper, large cavities are flooded with liquid and higher 
wall superheats are therefore required to nucleate the smaller 


cavities. 
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You et al [Ref. 7] investigated the onset of nucleate 
boiling of R-113 under atmospheric conditions using an 
electrically heated 0.13 mm diameter chromel wire and a 0.51 
mm platinum thin-film heater. Results of their repeated runs 
showed that the onset of nucleate boiling was very unstable 
and a probabilistic plot was utilized to display their 
results. They attributed the unstable behavior to the 
excellent wetting characteristics of R-113 and the poor 
reproducibility of contact angle. 

3. Contact Angle 

Contact angle between a liquid and a solid surface is 
one of the most important factors in phase change heat 
transfer. Griffith and Wallis [Ref. 27) concluded that 
contact angle is important in bubble nucleation primarily due 
to its effect on the stability of the bubble within the 
cavity. This importance was restated by Tong et al. [Ref. 26] 
and Cornwell [Ref. 28], since contact angle characterizes the 
wettability of a certain liquid on a specific surface. 

4. Surface Roughness 

Significant advances have been achieved in the 
development of enhanced boiling surfaces during the past 20 
years. Some of the commonly-used enhanced surfaces are shown 
in Figure 2.6 from Thome [Ref. 20]. Commercially available 
porous coatings and grooved/finned surfaces are capable of 


providing significant increases in heat transfer for various 
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surface/fluid combinations. Bar-Cohen [Ref. 9] states that 
porous coating provide a significant increase in the number of 
active nucleation sites per unit area and shift the active 
nucleation site population toward a larger mean radii by 
forming reentrant cavities, thus resulting in a lower 
incipient wall superheat. 

Numerous investigations have been conducted on 
enhanced surfaces and some are presented here. Chowdhury and 
Winterton [Ref. 29] report that an increased number of 
nucleation sites results in lower wall superheats, moving the 
nucleate boiling curve toward the left. Venart et al. [Ref. 
30) did a comparison using R-11 between a non-porous and two 
porous surfaces and found that the porous surfaces showed a 
significant decrease in wall superheat. Memory and Marto 
{[Ref. 8], and Wanniarachi et al. [Ref. 31 & 32] from their 
investigations using R-114 on smooth and enhanced tubes, 
concluded that the onset of nucleate boiling for an enhanced 
tube occurred at lower heat fluxes and wall superheats when 
compared to a smooth tube. 

5. Subcooling 

Surface subcooling is defined by Bergles and Chyu 
[Ref. 13] as the temperature difference between saturation 
temperature at the surface and the temperature of the bulk 
pool before the test. They found that the onset of nucleate 


boiling wall superheat for a smooth surface in R-113 increased 
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Pigure 2.6 Commercially available enhanced surfaces 
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with an increase in subcooling. Marto and Lepere [Ref. 14] 
tested three different tubes using different surface 
preparations and found that the ONB wall superheat increased 
with increasing subcooling. Furthermore, they found that pre- 
boiling followed by 30 minutes of subcooling has almost the 
same effect on the incipient wall superheat as subcooling the 
specimen overnight. Faw et al. [Ref. 33], using heating 
elements made of platinum wire, investigated the effects of 
pressure on the initiation of subcooled pool boiling in water. 
They used saturation temperatures of 100, 150, and 195 “Cc and 
subcoolings as much as 30 K. They observed that with 
increased subcooling, the ONB wall superheat also increased. 
This behavior was confirmed recently by Normington et al. 
(Ref. 34], using D80 and HT110. All of the above findings can 
be attributed to condensation of the vapor trapped in the 
cavities when subcooled, thus eliminating or reducing many 
potential active nucleation sites. 
6. Pressure and Temperature 

Pressure and temperature both play vital roles in the 
incipience of nucleate pool boiling. In 1964, Bergles and 
Roshenow [Ref. 35] developed the following equation for 
incipient wall superheats for liquids operating far from their 


thermodynamic critical point: 


(T, ~ Trae), = SloT, 3] (2.5) 
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They proposed using the bracketed fluid property group as a 
nucleation parameter. As pointed out by Bar-Cohen [Ref. 9] 

it can be seen that at saturation conditions, increasing 
system pressure and temperature will decrease the value of the 
above property group. This is not readily obvicus, but 1f we 
take R-114 for example, as saturation temperature increases 


the latent heat of vaporizatior (h,,), the liquid-vapor 


y 
specific volume (v,,) and the surface tension all decrease, 
thereby decreasing the wall superheat. Just like all other 
parameters affecting incipier se, several studies have been 
conducted on pressure-temperature effects. Eddington and 
Kenning [Ref. 36] found from their study that applying excess 
pressures of 1/3 bar to a system at 1 bar reduces the number 
of nucleation sites of a given size by about 50%. 

Hooper and Abdelmessih [Ref. 37] investigated the 
effect of decompression on nucleate boiling and noted that 
boiling initiation took place at higher superheats. This 
finding was confirmed by Weisman et al. [Ref. 38], Lienhard et 
al. [Ref. 39], and Faw et al. [Ref. 33]. 

7. Oil addition 

In most applications involving chillers, the mixing of 
the refrigerant with a miscible refrigeration oil is 
inevitable due to the use of hermetically sealed centrifugal 
compressors. Oil seals deteriorate with time and oil leaks 


into and mixes with the working fluid. Several studies have 
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been done on the effect of oil on nucleate pool boiling of 
refrigerant-oil mixtures. Stephan [Ref. 40] was probably one 
of the earliest investigators. He studied the effect of oil 
on the boiling heat transfer coefficient using R-12 and R-22 
on a smooth horizontal plate. He reported up to a 50% 
reduction in the boiling coefficient with 9% oil by mass, 
while 50% oil caused a 90% reduction. Daugherty and Sauer 
{Ref. 41] investigated the effects of R-11/oil and R-113/oil 
mixtures on boiling hysteresis from smooth tubes and found 
that o11 additior (up to 10%) caused a significant increase in 
the incipient wall superheat for both fluids under the same 
operating conditions. Sauer et al. [Ref. 42] investigated the 
effect of oil (up to 10%) carryover on the boiling performance 
of R-11 with commercial copper 19 fpi finned tube. They found 
that incipient wall superheat increased with increasing oil 
concentration and that the effect of oil on finned tubes was 
similar to smooth tubes. Wanniarachi et al. [Ref. 32] studied 
the effect of oil on the performance of R-114 from smooth and 
porous surfaces and concluded that the presence of oil 
Significantly delayed the ONB for the porous coated tube, but 
only had a small effect for the smooth tube. Memory and Marto 
{Ref. 8) in their study of R-114/oil mixtures, found that 
oil concentration had no systematic effect on the incipient 
point for smooth and low integral-fin tubes, but delayed it 
for structured type surfaces. Additionally, they concluded 


that boiling incipience is highly unpredictable and can easily 
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be triggered by an external source. Bertsch [Ref. 6] during 
his study of R-124/oil mixtures where he conducted many 
repeated experiments, concluded that oil addition does delay 
the point of incipience for a smooth tube. However, there was 
little measurable effect for a HIGH FLUX tube. This 
contradictory evidence is probably due to the unpredictable 


and statistical nature of incipience. 
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III. EXPERIMENTAL APPARATUS 


A. SYSTEM DESCRIPTION 

A schematic diagram representing the single-tube apparatus 
is shown in Figure 3.1, and a photograph is shown in Figure 
3.2. Detailed discussions on its design, construction and 
assembly were done by Karasabun [Ref.43]), and Reilley [Ref. 
44}. The apparatus was relocated to its existing location and 
calibrated by Sugiyama [Ref. 45]. Existing Pyrex glass Tees 
for the condenser and evaporator were replaced by Bertsch 
{[Ref. 6] with fiberglass re-inforced Pyrex glass Tees to 
withstand the higher operating pressure of R-124. 

The present system consists of the following components: 
(1) two identical re-inforced Pyrex-glass Tees, one used as 
the evaporator to boil the R-124, and the other used as the 
condenser to condense the refrigerant vapor; (2) a liquid 
refrigerant reservoir; (3) a refrigerant/oil subsystem; (4) an 
aluminum/plexiglass framework; (5) an ethylene-glycol/water 
mixture cooling system; (6) a vacuum pump; and (7) a data 
acquisition and instrumentation system. 

R-124 was boiled in the evaporator and was condensed in 
the condenser on a copper coil cooled by a refrigerated water 
and ethylene glycol mixture. Condensation rate was controlled 


via valve VC (see Figure 3.1). The condensate then returned 
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Figure 3.1 
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Schematic Diagram of Pool Boiling Apparatus 
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to the evaporator by gravity. A copper deflector (see Fig. 
3.3) was fitted above the condensate discharge tube at the 
bottom of the evaporator to prevent any vapor bubbles that 
might come from the discharge tube from interfering with the 


boiling tube. 


B. BOILING TEST SECTION 
1. Evaporator 
The evaporator was a Corning Pyrex fiberglass- 
reinforced glass Tees section designed to withstand a working 
pressure of 50 psig. Bertsch [Ref. 6] hydrostatically tested 
the system to 75 psig prior to the initial charging of R-124. 
Two aluminum endplates, 210 mm in diameter and 12.7 mm thick 
with thermocouple wells, oil entry fittings and Teflon inserts 
were bolted to the evaporator by means of cast-iron flanges 
and gasket assembly as shown in Figure 3.4 and Figure 3.5. 
2. Evaporator tubes 
Each evaporator tube tested was secured in place by 
two Teflon inserts and sealed by O-rings in the evaporator 
endplates. Strong-backs as shown in Figure 3.6 prevented the 
Teflon inserts from popping out due to internal pressure of 
the system. A schematic drawing of an evaporator tube is 
shown in Figure 3.7. 


Each evaporator tube was 450 mm in length and 15.9 mm 


in diameter. The middle portion of the tube was heated by a 











Figure 3.3 Photograph of Copper Deflector Plate 
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1 kW, 240-Volt stainless-steel-sheathed cartridge heater, 
6.35 mm in outside diameter, 203.2 mm in actual length with a 
heated length of only 190 mm. Memory and Marto [Ref. 8], and 
Bertsch (Ref. 6] provide a detailed fabrication process of an 
instrumented tube. 

Four tubes were tested in this study: (1) smooth tube, 
(2) GEWA-K 19 fin per inch (fpi) tube, (3) GEWA-K 26 fin per 
inch tube and (4) HIGH-FLUX porous tube. All the tubes were 
made of copper except the HIGH-FLUX tube that had copper 


particles on a copper-nickel substrate. 


C. CONDENSER SECTION 

The condenser was identical in size, shape and material to 
the evaporator shown in Figure 3.4. The refrigerant was 
condensed on a 3/8 inch outside diameter helical copper tube 
that was coiled inside the condenser. The condenser was 
mounted vertically as shown in Figure 3.1 with top and bottom 
plates bolted to the fitted cast-iron and gasket assembly. 
The top plate has both coolant entry and vacuum/charging 
fittings; the bottom plate was connected to the evaporator by 
a gravity feed return line via valve V5. The condenser was 
connected to the vapor outlet of the evaporator by a 2-inch 


diameter L-shaped aluminum piping. A relief valve (set at 50 


psig) and a Bourdon pressure gage were mounted on the tube. 








a) Corning Pyrex Glass Evaporator (D x f= %02x51 mm, 
L = 178 mm, = 127 mm) 





b) Cast Iron Flange and Gasket (d,= 190mm, d y= 210mm, 
L,= lumm, A = 21") 


Figure 3.4 Sketch of Pyrex Glass Vessels and Flanges 
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Figure 3.6 Endplate strongback to secure Teflon inserts 
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D. OIL ADDITION SECTION 

The oil addition section consisted of a 152.4 mm diameter 
by 152.4 mm high oil reservoir, and a 25.4 mm diameter by 355 
mm high graduated glass cylinder mounted above the evaporator 
as shown in Figure 3.1. Replenishing oil in the reservoir was 
achieved by removing the filler cap, and pouring oil through 
a funnel. The graduated cylinder was connected to the oil 
reservoir through valves V4 and V2. A miscible alkylbenzene 
oil, Zerol 300, was used with R-124. Oil concentrations (by 
weight) of 0, 3, and 10% were used for each of the tubes. 
Addition of oil into the refrigerant was accomplished with the 
following procedure: 

1. Ensure pneumatic hose is connected to the reservoir’s 
filler tube and charged. Flow of oil into the graduated 
cylinder from the cylinder is controlled by means of 
valve V2. 

2. Ensure refrigeration plants are on line and sump 
temperature is at least -10 °C. Start the condenser pump 
and lower the system pressure to about 5 psig by 
controlling condensate flow by means of valve VC. 

3. Open valve V4 to equalize the pressure in the graduated 
cylinder and in the evaporator/condenser (indicated by no 
more pressure change with the installed pressure gage). 
Secure the condenser pump. 

4. Pressurize the entire system by applying heat flux to the 
refrigerant pool by aligning the variac ring to 
approximately 90 KW. Secure heat flux when system 


pressure reaches approximately 30 psig. Then close valve 
v4 to contain the pressure in the tube. 


4. Start condenser pump and open valve vc fully to lower 
pressure in the system to approximately 5 psig. Then open 
valve V1 to allow the measured oil in the graduated 
cylinder to flow into the pool. 











Tne above procedure will expedite the o:1 addition 


process, the higher the pressure difference between the trop 


and bottom of the cylinder, the quicker the flow of o1. 


E. COOLING SYSTEM 
1. Water/Ethylene-Glycol Coolant Sump 
A plexiglass tank, 1.3 cm thick, 47.85 cmon its sides 
and 66.14 cm high contained a mixture of 13 gallons of 
ethylene-glycol and 25 gallons of distilled water. It was 
mounted at the base of the frame on a wooden platform, and was 
insulated with 22 mm thick rubber sheet. Supply and discharge 
ports were connected for the condenser and sump pumps. A 
thermocouple to measure sump temperature was irstalled through 
the condenser return line located on the top of the tank. 
2. Refrigeration Plants 
A 1/2 ton R-502 and 1/4 ton R-12 refrigeration system 


cooled the water/ethylene glycol mixture in the sump. Both 


systems consisted of the following components: (1) 
hermetically sealed compressor assembly, (2) air-cooled 
condenser, (3) receiver, (4) filter-dryer, (5) pressure 
regulator, (6) temperature control switch, and (7) 
thermostatic expansion valve. Each refrigeraticn plant was 


controlled by a temperature control switch and thermostatic 
expansion valve. The R-12 plant used a 9.5 mm coiled copper 


tubing installed in the sump tank as an evaporator as shown in 


Figure 3.1, while the R-502 plant used a counter-flow heat 











exchanger as an evaporator as shown in Figure 3.8. The sump 
temperature was maint iined at approximately -17 °C when both 
R-12 and R-502 plants were used. Figure 3.9 is a photograph 
of the R-12 and R-502 refrigeration plants. 
3. Pump and Control Valve 

A positive displacement, 8 gpm, 115V Burks turbine 
type pump was used to pump the ethylene glycol/water coolant 
mixture through cross-flow heat exchanger of the R-502 
refrigeration plant. A similar pump was used to pump the 
coolant mixture to the condenser via valve VC as shown in 
Figure 3.1. Valve VC was used by the operator to control the 
refrigerant pool saturation temperature and condensation rate. 
A discharge bypass line with valve V9 was installed from the 
pump to prevent overloading the positive displacement pump 
during low heat flux application (i.e. low coolant flow 


through the condenser). 


F. REFRIGERANT RESERVOIR 

To facilitate tube removal without dumping the refrigerant 
into the atmosphere, an aluminum reservoir 228.6 mm in 
diameter by 254 mm high was mounted just above the evaporator 
{see Figure 3.1). The reservoir was equipped with a sight 
glass for liquid-level monitoring. The reservoir can be 
connected or isolated to or from the system by means of valves 


V6, V7, Vl2 and V13. 
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G. FRAME 

An aluminum frame 107 cm X 51 cm X 61 cm contained all 
the components of the system with the exception of the R-502 
and R-12 refrigeration plants. The top portion of the frame 
contained the evaporator, concenser, oil addition system, and 
the refrigerant reservoir. For safety considerations and to 
help maintain the evaporator and condenser in an isothermal 
condition, all sides of the frame were covered with 12.7 mm 
plexiglass. The left and right side covers were hinged to the 
frame to provide easy access to the assembly. The bottom and 
top portion of the frame were covered with an aluminum plate 
and plywood respectively. The coolant sump tank was mounted 


at the base of the frame. 


H. INSTRUMENTATION 
1. Power Measurement 
A 240 volt AC source, adjustable from 0 to 220 volts 
and 0-5 amps by a variac controller, was used as the power 
source for the heater in the boiling tube. The AC current and 
voltage RMS sensors converted the AC input into DC voltages 
that were then used as input to the Data Acquisition and 
Control Unit as shown in Figure 3.10. 
2. Temperature Measurement 
Temperatures were measured using the thermal emf from 
30 gage (0.25 mm diameter) copper-constantan (type T) 


thermocouples. Eight thermocouples were spaced both 
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Figure 3.9 Photograph of refrigeration plants 
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circumferentially and axially over the heater sleeve, as shown 
in Figure 3.11, to measure the tube wall temperatures. Three 
thermocouples inserted into the evaporator thermocouple wells 
measured the evaporator pool temperatures. A single 
thermocouple inserted into the coolant sump tank measured the 


sump temperature. 


I. DATA ACQUISITION AND REDUCTION PROGRAM 

Bertsch’s [Ref. 6] DRPGB data acquisition and reduction 
program was used in this study. The interactive data 
acquisition program directed the Hewlett-Packard 9826 computer 
and a Hewlett-Packard 3497a Data Acquisition Unit to read and 
process all thermocouples, current and R.M.S. voltage 
readings. Raw data were stored in the computer’s hard drive, 
and a hardcopy printed on a Hewlett-Packard Inkjet printer. 


Channel assignments are listed in Table 3.1. 


J. PROCEDURE FROM KEYBOARD 
The following interactive procedure was used for data 
acquisition: 


1. Load data acquisition/reprocessing program by 
pressing F5 (load button) then type in ‘DRPGB’. 


Press F3 (run button) to run the program. Select 
‘0’ (TAKING DATA OR RE-PROCESSING PREVIOUS DATA 
OPTION) . 


2. After the prompt "GIVE A NAME FOR THE RAW DATA 
FILE", type filename of data to be processed and 
stored (see Appendix A for data file naming scheme 
used). 
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3. Enter any defective thermocouple locations that were 
previously identified by program SETUP71 (see 
Appendix © of Bertsch [Ref.6]). 


4. Select test tube after the prompt "SELECT TUBE 
NUMBER" by typing the number corresponding tube 
type displayed on the screen. For example, type ‘5’ 
for HIGH-FLUX tube, ‘7’ for GEWA-K 26 fpi, ‘13’ for 
GEWA-K 19 fpi. 


5. Select option ‘1’ (SET DESIRED HEAT FLUX), then 
adjust variac to the desired heat flux value. 


6. Adjust coolant flow through the condenser with valve 
vc as required to obtain the desired pool and vapor 
temperatures (TL1, TL2, and TLV) close to saturation 
temperature (2.22 + .01 CC). 


7. When saturation temperature was reached, and after 
waiting 5 minutes for the temperature to stabilize, 
press Fl to display the select option, then press 
‘0’ to take data. This will direct the date 
acquisition computer to scan all the channels as 
listed in Table 3.1, do the required calculations, 
display them on the screen and print the results. 


8. After the result is printed, press ‘1’ after the 
prompt "OK TO STORE THIS DATA SET" to store data 
into the file name created in step 2. 


9. -After the prompt “WILL THERE BE ANOTHER’ RUN 
(1=Y,0=N)?", press ‘1°. 


10. Repeat steps 5 to 9 for each change in heat flux. 


For the last heat flux entered, press ‘0’ for step 
no. 9. 
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Bottom 


(a) View of the boiling tube thermocouple locations as seen 
from the front of the experimental apparatus. 


Bottom 
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(b) Thermocouple sleeve unwrapped (act section A-A) to show 
the relative locations of the thermocouple channels 
(aIl dimensions in millimecers). 


LEFT-END VIEW RIGHT-END VIEW 
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(c) End views of the boiling tube. 


Figure 3.11 Thermocouple locations of instrumented tube 
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Iv. EXPERIMENTAL PROCEDURE 


A. SYSTEM PREPARATION 
1. Vacuum Test 
Once the evaporator tube was installed, and prior to 
charging with refrigerant, the system was evacuated by opening 
valves V1l1 and V8 and activating the vacuum pump (see Figure 
3.1) to approximately 29 inch Hg. If no leak was detected 
within 30 minutes (i.e. no drop in vacuum), the system was 
charged with R-124. 
2. Pressure Test 
If there was a significant vacuum loss during the 
vacuum test, the system was pressurized with air to 20 psig 
via valves V10 and V8. Leaks were detected using soap/water 
solution. After identifying and correcting the leak source, 
the system was again evacuated and step 1 was repeated. 
3. System Charging with R-124 
After successful system integrity was accomplished, 
the evaporator was charged with refrigerant to a level of 20 
mm above the top of the test tube using the following 
procedure: 
a. Gravity Fill from the Reservoir 
Valves V12 and V7 were opened to equalize pressure 


between the reservoir and the evaporator. Then valves V6 and 


47 


V13 were opened to allow the liquid refrigerant to drain into 
the evaporator. At the required liquid level in the 
evaporator, valves V6 and V7 were closed to isolate the 
reservoir. 
b. From Refrigerant Cylinder 

The cylinder was first connected to valve V10 by 
means of a charging hose. Valves V5, V8 and V10 were opened 
and refrigerant was released from the cylinder by opening the 
liquid valve. Slightly before the required liquid level was 
attained, the cylinder charging valve was closed, followed by 
the closure of valves V10 and V8. 

4. Degassing Procedure 
a. Degassing Procedure A 

After charging the system with refrigerant, a leak 
test waS conducted using the hand-held “Halogen Leak 
Detector", by energizing the detector and placing its probe 
near all flange connections and hose fittings. (If a leak is 
detected, the detector will emit a series of beeping sounds.) 
After correcting the leak (i.e. tightening flange bolts or 
hose fittings), the system was degassed by first boiling the 
refrigerant at a power setting of approximately 90 KW (on the 
variac ring) for 10 minutes with the condenser pump running to 
prevent overpressurizing the system. Any noncondensable gases 
produced by this process were collected in the top portion of 


the condenser. The heat flux and condenser pump were then 
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secured and the system allowed to sit for approximately 10-20 
Minutes; this allows the condensed refrigerant to drain back 
into the evaporator, leaving just the noncondensable gases. 
Finally, the vacuum pump was turned on and valves V1l1 and V8 
were opened for approximately 10 seconds to suck out the non- 
condensable gases. 
b. Degassing Procedure B 
This degassing procedure came about as a result of 
consistently low incipient superheats after long settling time 
(time between runs of approximately 12 hours or more 
particularly on the HIGH FLUX tube). This procedure was very 
similar to degassing procedure A except for the following 
differences: (1) the pool was boiled for 30 minutes instead of 
10 minutes, (2) the system was secured for approximately 30 
minutes prior to activating the vacuum pump, (3) the vacuum 
was turned on for 20-30 seconds. Procedure B was also 
conducted after every 011 addition due to the visible presence 
of entrapped air in the oil in the graduated cylinder prior to 
oil addition. 
5. Data Acquisition Channel Check 
After connecting all evaporator tube thermocouples, 
and inserting the three pool thermocouples into their 
respective wells, program SETUP71 (see Bertsch (Ref. 6]) was 
loaded and run to check the output of all thermocouples. Any 


inconsistent readings from the tube thermocouples were then 
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noted and excluded from the data acquisition portion. For any 


other erroneous readings, appropriate action was taken. In 
some cases, the apparatus was secured for at least 6 hours 
(usually overnight) to allow the tube surface to become fully 


wetted. 


B. NORMAL OPERATING PROCEDURE 
The following procedure was utilized to conduct test. 
This was closely followed from Bertsch [Ref. 6]: 


1. Both R-502 and R-12 refrigeration systems were operated 
for approximately 2 hours to cool the sump to -17 °C. 


2. The data acquisition unit, computer and variac control 
panel were energized. 


3. Program SETUP71 was loaded and run. All data acquisition 
channels were checked, and the refrigerant pool 
temperature was slowly reduced to 2.2 °C by circulating 
coolant through the condenser via valve VC. 


4. Data acquisition program DRPGB was loaded and run. 
5S. An initial heat flux (600-700 W/m’) was entered in the 


program and the variac was adjusced to this preset value. 


6. Coolant flow through the condenser was adjusted by 
opening slightly valve vV10 until desired saturation 
temperature of 2.22 + 0.1 °C was achievea. 


7. Data were collected after the system had stabilized for 
5 minutes. 


8. Heat flux was increased logarithmically (i.e. in equal 


sized increments on a logarithmic plot) and steps 6 and 
7 were repeated. 
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V. DISCUSSION OF RESULTS 


A. REPRODUCIBILITY 

Reproducibility/repeatability runs were initially 
conducted for familiarization with the equipment = and 
operating/testing procedures. Figures 5.1 and 5.2 are plots 
of increasing heat flux runs for the two degassing procedures 
taken on four different days. Good agreement tk:oughout, 
including the onset of nucleate boiling (indicated by arrows) 
is evident. Small amounts of scatter at lower heat fluxes 
(from 500 to 3000 W/m?) which disappears at mid-range heat 
fluxes (3000 to 6000 W/m’) can be seen in both figures. This 
scatter is probably due to thermal gradients within the pool 
in the convection region where the mixing is not very 
efficient. Large local variations in the wall and pool 
thermocouple readings tend to support this. For example, the 
run conducted on 29 Jan 93, at an initial heat flux of 430 
W/m?, three of six wall thermocouples (located close to both 
ends of the heated length of the tube) and one of the three 
pool thermocouples read about 5% lower than the other 
thermocouples. The same thermocouples read 2% lower at a heat 
flux of 1750 W/m’. In the nucleate boiling regime (6000 to 
90,000 W/m*), excellent repeatability was achieved at all heat 


fluxes. The corresponding decreasing heat flux runs, shown in 


St 








Figure 5.3 for degas procedure A also displays a very 
repeatable behavior. In the boiling region, the pool and wall 
temperature are much more steady due to better mixing in the 
pool. 

Figure 5.4 is a comparison of a representative increasing 
and decreasing heat flux run with the data of Bertsch [Ref. 6] 
for the smooth tube. As before, curves #1 (present data) and 
#3 (Bertsch) show slight scatter in the low heat flux natural 
convection region but merge at medium heat’ flux. Curve #1 
nucleates at a wall superheat of 15.4 °C, while curve #3 
nucleates at a lower wall superheat of 8.5 °C. The difference 
in incipient wall superheat results ina noticeable difference 
in the starting point of the nucleate boiling region. Curve 
#1 has a thermal overshoot (thermal overshoot is taken to be 
the horizontal jump between the incipient point and the 
decreasing heat flux curve (Memory and Marto [Ref. 8]) of 
about 11 °C while curve # 3 yielded only a 4 °C thermal 
overshoot. The reason for this difference is not clear since 
the tests were, in essence, identical. At higher heat fluxes 
in the nucleate boiling region, good agreement was obtained. 
In the decreasing heat flux runs, curves #2 and #4 were almost 
identical. 

Other comparisons and repeatabilities are shown in Figures 
5.5 to 5.15 for various tubes and various oil concentrations. 
In Figure 5.5, with 3% oil, the present data show a much 


higher beat transfer coefficient (lower wall superheat 
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throughout the natural convection region) and a lower ONB 
compared to Bertsch [Ref. 6]. Good repeatability was again 
shown in the nucleate boiling region. In the decreasing runs, 
both data are very repeatable at heat fluxes between 90,000 
W/m to 30,000 W/m*?; but separation occurred at lower heat 
fluxes. As mentioned earlier, reasons for the differences are 
not clear. Figure 5.6 with 10% oil, show the same trend in 
the natural convection and nucleate boiling regions between 
Bertsch [Ref. 6] and the present data, except that a 
relatively higher ONB was obtained. The decreasing runs 
display good level of repeatability. 

Figures 5.7, 5.8, and 5.9 are GEWA-K 19 fpi finned tube 
comparison plots of pure R-124, 3% oil and 10% oil data of 
Bertsch [Ref. 6] and the present study respectively. All 
plots display good repeatabilities throughout the natural 
convection and nucleate boiling regions, with consistently 
higher ONB points obtained from the present study. This may 
be attributed to the difference in degassing procedures used. 

Figures 5.10, 5.11, and 5.12 are GEWA-K 26 fpi finned tube 
comparison plots for pure, 3% and 10% oil respectively of 
Bertsch [Ref. 6] and the present study. All plots show 
Significant differences throughout in the natural convection 
region and most notably in the nucleate boiling region where 
a consistent difference of about 10 K can be seen at 100,000 
W/m’. This can probably be due to a better cleaning procedure 


used in the present study (cleaning procedure B as discussed 


ee) 











in Chap IV) in which the fin gaps were brushed with alcohol 
then with acetone before installation, vice just spraying with 
acetone (cleaning procedure A) used by Bertsch [Ref. 6]. 
Figures 5.13, 5.14 and 5.15 are HIGH FLUX tube comparison 
plots. Consistently higher ONR points from the present study 
(due to better degassing procedure) and good agreement in the 


nucleate boiling region are evident. 


B. EFFECT OF OIL ON BOILING INCIPIENCE 

The effect of oil on the onset of nucleate boiling has 
been studied quite extensively in the past by Memory and Marto 
{Ref. 8], Chongrungreong et. al [Ref. 18], Stephan and 
Mitrovic [Ref. 22], Wanniarachi et. al. [Refs. 31 and 32], 
Stephan [Ref. 40], Daugherty and Sauer [Ref. 41], Sauer et al. 
(Ref. 42}, Karasabun [Ref. 43], Reilley [Ref. 44], Sugiyama 
(Ref. 45], and Bertsch [Ref. 6]. The effect of oil on 


incipient wall superheat can best be explained by the use of 


equation (2.5) and the corresponding nucleation parameter 


given by: 


(5.1) 


It can be seen from equation (2.5) that the incipient wall 
superheat depends on the parametric group inside the bracket 
as discussed by Bergles and Roshenow [Ref. 35]. Saturation 


temperature, surface tension and viscosity increase with 











increasing oil concentration as pointed out by Bar-Cohen [Ref. 
9]. Sauer et al. [Ref. 42] said: 

Several thermodynamic properties of the refrigerant 

are significantly changed by the addition of oil. The 

vapor pressure of the refrigerant-o1il mixture is lower 

than that for the pure refrigerant at the same 

temperature. Both viscosity and surface tension are 

increased when oil goes into solution with the 

refrigerant....An increase in surface tension 

decreases bubble formation for a given wall superheat. 

Increased liquid viscosity tends to reduce turbulence 

and dampen out eddy formations in the liquid flow to 

the surface. These may be the principal factors which 

account for the various effects of oil on nucleate 

boiling. 
In the present study, however, there are no available 
thermodynamic properties of the oil (Zerol 300) to determine 
the mixture’s latent heat of vaporization, specific volume 
and other thermodynamic properties, therefore, it is 
impossible at this time to predict precisely the incipient 
boiling superheat with the use of equations (2.5). 

1. Smooth Tube 

Figures 5.16 5.17, and 5.18 are plots of seven 

increasing heat flux runs for pure R-124, R-124/3% oil mixture 
and R-124/10% oil mixture respectively. These figures 
demonstrate the scattering of the data in the natural 
convection region, particularly at lower heat fluxes due to 
reasons given earlier. In Figure 5.16, the data points tend 
to merge towards the Churchill-Chu [Ref. 46] correlation 


prior to the onset of nucleate boiling. The highest heat 


fluxes reached prior to boiling for 0%, 3% and 10% oil 
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concentrations were 12.2 Kw/m’, 12.2 Kw/m, and 14.9 Kw/m 


respectively. It should be noted that once nucleation 
occurred, very few data points in the nucleate region were 
taken as can be seen from the figures. For the few data 
points taken, excellent repeatability were always obtained. 
Apart from the repeatability/comparison data presented in 
Figures 5.1 to 5.15, no decreasing heat flux data were taken. 
For the pure R-124 case, incipient wall superheats ranged from 
30.2 K to 35.3 K with corresponding thermal overshoots from 
21.0 K to 25.8 K. For the 3% oil/R-124 mixture, incipient 
superheats were from 31.2 K to 37.3 K with thermal overshoots 
of 24.0 K to 28.7 K. The 10% oi1/R-124 mixture resulted in 
incipient superheats of 36.2 K to 40.1 K with thermal 
overshoots ranging from 27.0 K to 30.1 K (see Table 5.1). 
Hence, ONB was fairly consistent and only varied typically by 
5 °C for numerous tests runs under identical conditions. In 
Figures 5.17 and 5.18, the Churchill and Chu correlation is 
not given as it is strictly valid for pure substances. 
However, the data still tend to converge as heat flux is 
increased in the convection region for similar reasons as 
those already given. 

Figure 5.19 is a probabilistic summary of the onset of 
nucleate boiling. This form of presentation for incipient 
superheat was introduced by You et al. [Ref. 7] as the most 
convenient way to display the randomness behavior of the onset 


of nucleate boiling. Figure 5.19 shows the incipient wall 


56 








(a) ebuey 
BAOOUSIBVAQ [RuUrouL 


XxN1Id HDIH - JH 
TdJ 92 M-WMED - 9% ID 
TdJ 6T N-WMAD - 6T AD 
HLOOWS - WS ‘S.LON 


(a) ebuey | (,0/Mx) eoueTdzouy 
qQeeyzedns QueTdtfour zofy1zd peutejqqAe 





NOSTAYWdWOD ONITIO€ ALVSTIONN JO LESNO T°S STQRL 





57 





superheats attained for the 7 tests conducted with pure R-124, 
3% and 10% oil concentrations as percent probability of 
nucleation occurring. For example, for the 3% oil curve, all 
seven runs reached 31.5 °C without nucleating. Three of the 
seven runs nucleated at a wall superheat less than 32.8 °C 
and four of seven reached nucleation point at a higher wall 
superheat; this means that the probability of nucleation 
occurring at 32.8 °C is between 3/7 and 4/7 or 3.5/7 (50%). 
At a wall superheat of 37.3 °C all seven runs reached the 
nucleation point so that a probability of 100% exists. 
Notice that with this probabilistic representation, it is easy 
to see that the presence of oil delays the onset of 
nucleation. Figure 5.20 is a comparison between the smooth 
tube data of Bertsch [Ref. 6] with the present data. The 
difference can be attributed to factors such as degassing 
(which is discussed in more detail later) and data gathering 
techniques, by the two investigators. One thing in common, 
however, is that increasing 0il concentration definitely seems 
to result in an increase in incipient wall superheat, but 
the effect may be small in comparison to the large 
probabilistic variations due to operating procedures. 
2. GEWA-K 19 Tube 

Figures 5.21, 5.22, and 5.23 show the effect of oil 

for seven increasing heat flux runs for pure R-124, R-124/3% 


oil mixture and R-124/10% oil mixture respectively while 
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boiling from a 19 fpi finned tube (GEWA-K). As before, data 
scattering was evident at low heat fluxes, with a merging 
trend at mid-range heat fluxes. The Churchill/Chu correlation 
is also given in Figure 5.21 and appears to be well to the 
right of the data. However, the reason for this is due to the 
choice of tube surface area in the heat flux calculation, 
where root diameter has been used. Memory and Marto [Ref. 8] 
have shown that if the actual wetted area were used in the 
calculation of heat flux, then the data would drop back very 
close to the correlation. Therefore, one must be careful when 
interpreting such results. The maximum heat fluxes attained 
prior to ONB were 12.3 kW/m’, 22.5 kW/m* and 24.5 kW/m’ for 0%, 
3% and 10% oil respectively. For pure R-124, incipient wall 
superheats ranged from a minimum value of 9.4 K to a maximum 
of 11.6 K with thermal overshoots ranging from 7.2 K to 9.4 K. 
These are significantly lower than the smooth tube, 
demonstrating the improved heat transfer performance from 
these enhanced tubes. At 3% oil, the incipient wall 
superheats were more dispersed from 10.0 K to 19.6 K with 
thermal overshoots ranging from 7.3 K to 16.4 K. Increasing 
the oil concentration to 10% yielded incipient wall superheats 
ranging from 14.0 K to 20.2 K with thermal overshoots from 
11.9 K to 16.3 K. Ona probabilistic basis, Figure 5.24 shows 
a wide variation in the incipient wall superheats for the 3% 


and 10% oil mixtures which re-affirms the randomness of ONB. 
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Also, the figure demonstrates the increase in ONB as oil 


concentration increases, as from the smooth tube. 
3. GEWA-K 26 Tube 

Figures 5.25, 5.26, and 5.27 show the results of oil 
on seven increasing heat flux runs for pure R-124, R-124/3% 
oil mixture and R-124/10% oil mixture respectively, while 
boiling from a 26 fpi finned tube (GEWA-K). Just like the 
smooth and 19 fpi GEWA-K plots, scattering at low heat fluxes 
is clearly evident. The maximum heat fluxes attained prior 
to ONB were 33.9 kW/m’, 28.6 kW/m? and 38.4 kW/m’ for 0%, 3% 
and 10% oil respectively. Incipient wall superheats were from 
17.5 K to 22.4 K with thermal overshoots ranging from 13.0 K 
to 17.3 K for pure R-124. These values were much larger than 
those noted earlier for the GEWA-K 19 fpi tube (see Table 
5.1). The reason for this difference is not known but may be 
due either to different surface roughness effects or different 
micro-convection effects within the channels between the fins. 
With 3% o11, the incipient wall superheats ranged from 12.9 K 
to 21.1 K, giving thermal overshoots of 9.2 K to 17.0 K. At 
10% oil concentrations, incipient wall superheats ranged from 
20.4 K to 23.7 K with corresponding thermal overshoots from 
14.6 K to 17.8 K. Figure 5.28 summarizes probabilistically 
the ONB data for 0%, 3% and 10% 011 concentrations. Note that 
for 3% oil, the data are to the left of the data for 0% oil. 


This shift was not obtained with the GEWA-K 19 fpi tube (see 
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Figure 5.24). However, this shift is very repeatable as these 
curves are data taken from 7 identical tests. The reason for 
the shift 1s not clear but could be a combination of changes 
in the thermophysical properties (when oil is added) together 
with favorable conditions for bubble growth from the vapor 
nuclei trapped in the cavities. It may be that the 26 fpi 
tube, with its narrow gap width between fins, is conducive to 
bubble growth whereas the 19 fpi tube is not; there may be a 
critical bubble diameter which can grow on the 26 fpi tube 
(due to its narrow channels) but which condenses and collapses 
on the 19 fpi and smooth tubes. This of course is conjecture 
and requires further thought. Further addition of oil (up to 
10% oil) causes a delay due to the change in nucleation 
parameter properties as noted in equation (5.1). 
4. HIGH FLUX Tube 

Figures 5.29, 5.30, and 5.31 present the seven runs 
conducted for 0%, 3% and 10% oil concentrations respectively 
for the HIGH FLUX tube. Lower incipient points and 
Significantly larger thermal overshoots were obtained for all 
oil concentrations (see Table 5.1). Incipient wall superheats 
ranged from 9.2 K to 11.5 K with resulting thermal overshoots 
of 8.6 K to 11 K for 0% oil. The scatter is slightly more 
evident for this tube in the convection region but is still 
converging as heat flux increases. The data fall very close 


to the Churchill/Chu correlation, demonstrating that in this 
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region, a HIGH FLUX tube behaves in a similar way to a smooth 
tube. At 3% oil concentration, lower incipient wall 
Superheats were obtained than for 0% oil (like the GEWA-K fpi 
tube) from 6.5 K to 11.2 K with corresponding thermal 
overshoots ranging from 6.7 K to 10.6 K. The 10% oil 
concentration showed a large increase of incipient wall 
Superheat compared to pure R-124 values. Incipient points 
ranged from 13.8 K to 21.6 K, and thermal overshoots ranging 
from 8.6 K to 11.6 K. Figure 5.32 displays the variation in 
incipient wall superheats in a probabilistic way. Note again 
the decrease in incipient superheat at 3% oil concentration, 
Similar to that found for the GEWA-K 26 fpi tube. It may be 
that the nucleating cavities on these tubes are similar. 
Figure 5.33 compares the incipient wall superheats of Bertsch 
[Ref. 6] with the present data for the HIGH FLUX tube for all 
oil concentrations. It can be seen that the present data have 
consistently much higher wall superheats than those of 
Bertsch. The reason for this is thought to be due to the 


different degassing procedures, which is discussed later. 


C. EFFECT OF SURFACE ON BOILING INCIPIENCE 

Figures 5.34 shows a comparison of representative data for 
GEWA-K 19 fpi and a smooth tube. It can be seen that in the 
natural convection region the GEWA-F. 19 fpi tube has a higher 
heat transfer coefficient (lower wall superheat) than the 


smooth tube. As mentioned earlier, this is because the heat 
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flux was computed using a surface area based on the tube’s 
root diameter instead of the actual finned area. Dividing the 
heat flux by 2.66 (ratio of the actual finned area to the area 
using the root diameter) shifts the curve closer to the smooth 
tube (and the Churchill/Chu correlation). An enhancement still 
1s evident, indicating that the presence of the fins affects 
the local convection currents around the tube. In Figure 
5.35, using an area conversion ratio of approximately four for 
the GEWA-K 26 fpi tube, the resulting actual finned area heat 
flux curve also behaves more like a smooth tube in the natural 
convection region. This is similar to that reported by Memory 
and Marto [Ref. 8] for R-114. Perhaps fin spacing is an 
important variable in studying these convective effects. 
Figures 5.36, 5.37 and 5.38 are probabilistic comparisons 
of ONB from each enhanced tube (GEWA-K 19, GEWA-K 26 and HIGH 
FLUX) with the smooth tube. As shown in Figure 5.36, both the 
GEWA-K 19 fpi finned tube and the smooth tube display an 
increasing incipient point as oil increases. At 50% 
probability (using the smooth tube data as the basis of 
comparison), the onset of nucleate boiling decreased by 
about 60% for all oil concentrations in comparing the 19 FPI 
finned tube to a smooth tube. Figure 5.37 shows a similar 
decrease of approximately 40% for all oil concentrations for 
the 26 FPI finned tube. As shown in Figure 5.38, with the 
HIGH FLUX tube, the incipient wall superheat decreased by as 


much as 75% when compared to a smooth tube. This more 
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pronounced decrease can be attributed to the increase in 
nucleation sites on enhanced surfaces and their greater 
propensity to nucleate at lower wall superheat due to their 


larger relative size. 


D. EFFECT OF DISSOLVED GASES 

Dissolved gases in the refrigerant can significantly lower 
the point of incipience as described by Bar-Cohen [Ref. 9], 
Murphy and Bergles [Ref. 24] and Sugiyama [Ref. 44]. The 
effect of dissolved gases can be seen in Figures 5.16 and 
Figure 5.39. Degassing procedure A (as described in Chapter 
IV) was used initially for all tubes tested. With degassing 
procedure B, incipient wall superheats of about twice those 
obtained using degassing procedure A were obtained. Figure 
5.40 is a probabilistic plot of the effect of degassing 


procedures for a smooth tube with 0% 3% and 10% oil 


- °° 


concentrations. In Figure 5.41, the effect of dissolved gases 
over a long period of settling time (time between runs) can be 
seen. For overnight settling time (greater than 10 hours), 
about 73% lower incipient wall superheats were obtained 
compared to those obtained with just 0.5 hours or greater 
settling time. It is reasonable to assume that dissolved 
gases over long settling time will tend to be absorbed back 
into the refrigerant pool. Figure 5.41 clearly shows the 


trend of earlier incipience for longer times between runs. 
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E. EFFECT OF PRESSURE 

The desire to be able to conduct more runs per day during 
the early part of the research led to pressurization of the 
system prior to a run. Bergles and Roshenow [Ref. 35], Bar- 
Cohen (Ref. 9], and Hooper and Abdelmessih [Ref. 37] 
discovered that creasing pressure above saturation pressure 
at operation tends to delay the incipient superheat. 
Pressurizing the system prior to a run deactivates more of the 
nucleation sites by flooding them out (i.e. rewetting the 
cavities). 

Pressurization was achieved by pre-boiling the pool using 
maximum heat flux (around 90 kW/m’) without running the 
condenser until a system pressure of 30 psig was reached. The 
heater was then secured and the pool was allowed to settle for 
the desired period. Figures 5.42 and 5.43 display the effect 
of pressurization on 0% and 3% oil/R-124 mixture respectively. 
The number of hours displayed in both figures corresponds to 
the settling time of the pool. It can be seen from these 
figures that there is a noticeable significant delay in the 


incipient wall superheat when the system was pre-pressurized. 


F. EFFECT OF TUBE CLEANING 
During the early runs conducted with the HIGH FLUX tube 
surface, the tube was cleaned (after tests with oil) following 


the cleaning procedure of Bertsch [Ref. 6] by spraying it with 


acetone and allowing it to air-dry before installation 
















(Cleaning Procedure A). These tests resulted in higher 
incipient wall superheats similar to the smooth tube for pure 
R-124, 3% and 10% oil concentrations. It seems that the use 
of alkylbenzene oil tends to clog the reentrant cavities of 
the HIGH FLUX tube making the tube behave more like a smooth 
tube. These effects were not found for the R-114 tests and 
Maybe the mineral oil used there was easier to clean off the 
tubes. A new cleaning procedure (procedure B) was established 


which entailed brushing the porous surface with a soft bristle 








brush (i.e. toothbrush) using ethyl alcohol and then with 







acetone. The tube was then allowed to air-dry prior to 






installation and resulted in incipient wall superheats that 


were approximately 50% lower than those obtained using 







cleaning procedure A. Figure 5.44 shows the difference 


between the two cleaning procedures for 0% oil runs on a HIGH 






FLUX tube. Figure 5.45 is a probabilistic comparison between 





the cleaning procedures for 0%, 3%, and 10% oil 







concentrations. 





G. EFFECT OF SUBCOOLING 






Figure 5.46 and 5.47 demonstrate the effect of subcooling 


on a smooth tube with pure R-124 and 10% oil/R-124 mixture 






respectively. Subcooling was conducted by lowering the pool 





temperature below saturation by opening valve VC (see Figure 






3.1) fully, thus allowing maximum flow of coolant through the 





condenser. The subcooling value is the difference between the 








Saturation temperature of the pool (i.e. 2.22 °C for R-124) 
and the measured pool temperature that can be obtained from 
running either SETUP71 program or the data acquisition program 
(DRPGB). It can be seen from Figure 5.46 that subcooling the 
pool below the saturation temperature for as long as one hour 
resulted in an incipient wall superheat close or even higher 
than the value obtained when the pool was allowed to settle 
overnight. This is because by subcooling, many of the 
potential nucleation sites are deactivated by condensation, 
Similar to what happens by leaving the pool to settle for long 


periods (i.e. rewetting). 
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Figure S.1 Smooth Tube 0% oil degassing procedure B 
repeatability runs 


68 


















tg 
| : 
4 a 
a 10 E ais 
i nN INCIPIENT 
= a < SOILING 
> a / 
4 @, H 
S| 
> 
= 
zi 
— 
< 
= 
10° 
| + 29 JAN 93) | 
| x 30 JAN 93 
o 31 JAN 93 
» 01 FEB 93 
Churchil!-Chu 
402 _! 








10°! 10° 10! 102 
WALL SUPERHEAT (CK) 


: Figure 5.2 Smooth Tube 0% oil, degassing procedure A 
increasing heat flux repeatability runs 
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Pigure 5.3 Smooth Tube 0% oil, degassing procedure A 
decreasing heat flux repeatability runs 
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Figure 5.4 Comparison of smooth tube data with 0% oil 
(increasing and decreasing heat flux) 
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Pigure 5.5 Comparison of Smooth tube 3% oil data 
(increasing and decreasing heat flux) 
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Pigure 5.6 Comparison of smooth tube 10% o1i1 deta 
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Figure 5.7 GEWA-K 19 fpi 0% oil data comparison 
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Pigure 5.8 GEWA-K 19 fpi 3% oil data comparison 
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Figure 5.9 GEWA-K 19 fpi 10% oil data comparison 
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Figure 5.10 GEWA-K 26 fpi 0% oil data comparison 
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Figure 5.11 GEWA-K 26 fpi 3% oil data comparison 
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Figure 5.12 GEWA-K 26 fpi 10% oil data comparison 
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Pigure 5.14 HIGH FLUX 3% oil data comparison 
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Pigure 5.15 HIGH FLUX 10% oil data comparison 
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Pigure 5.16 Smooth tube 0% oil incipience tests 
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Figure 5.17 Smooth tube 3% oil incipience tests 
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Pigure 5.18 Smooth tube 10% oil incipience tests 
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Figure 5.21 GEWA-K 19 fpi O% oil incipience tests 
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Figure 5.22 GEWA-K 19 fpi 3% oil incipience tests 
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Figure 5.23 GEWA-K 19 fpi 10% oil incipience tests 
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Figure 5.25 GEWA-K 26 fpi 0% oil incipience tests 


92 


HEAT FLUX (W/M°2) 





10° 


« 
—_— 
ws 
2 f 
104 2? 
vad 
@ 
x 
Xtand 
xO 
-@ 
0 
a 
of 
3 
10 Ga 
em a* 
i RUN NR. DELTA T (K) 
oO 1 13.9 
x 2 12.9 
° 3 13.6 
‘ 4 16.2 
) 5 16.7 
+ 6 19.5 
a 7 21.1 
102 
10°! 10° 10! 102 


WALL SUPERHEAT (K) 


Figure 5.26 GEWA-K 26 fpi 3% oil incipience tests 
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Pigure 5.27 GEWA-K 26 fpi 10% oil incipience tests 
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Figure 5.29 HIGH FLUX 0% oil incipience tests 
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Pigure 5.30 HIG. .1°"” 3% oil incipience tests 
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Pigure 5.31 HIGH FLUX 10% oil incipience tests 
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Figure 5.34 GEWA-K 19 f£pi actual finned area and root 
diameter comparison with smooth tube at 0% Oil. 
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Figure 5.35 GEWA-K 26 fpi actual finned area and root 
diameter comparison with smooth tube at 0% oil 
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Pigure 5.39 Smooth Tube 0% oil, degassing procedure A 
incipience tests 
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Figure 5.41 Effect of dissolved gases on incipient wall 
superheats on HIGH FLUX tube (varying settling time) 
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Figure 5.42 Pressure effect on incipient wall superheats on 
a smooth tube with 0% oil 
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Figure 5.43 Pressure effect on incipient wall superheats on 
a smooth tube with 3% oil 
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Figure 5.44 Cleaning procedure effect comparison on 
incipient wall superheat on a HIGH FLUX tube with 0% oil 
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Pigure 5.46 Subcooling effect on incipient wall superheat 
from a smooth tube with pure R-124 
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Figure 5.47 Subcooling effect on incipient wall superheat 
from a smooth tube with 10% oil/R-124 mixture 
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VI. CONCLUSIONS 


Onset of nucleate boiling 1s very random. 


An oil concentration of 3% tends to delay the onset of 
nucleate boiling (compared to pure R-124) on a smooth 
tube and GEWA-K 19 fpi tube, but tends to lower the onset 
of nucleate boiling on a GEWA-K 26 fpi and HIGH FLUX 
tube. The reason for this is not precisely known, but is 
very repeatable. 


An oil concentration of 10% delayed the onset of nucleate 
boiling for all tubes. 


Presence of dissolved gases in the refrigerant tends to 
lower the onset of nucleate boiling by increasing the 
number of entrapped vapor nuclei of the surface. 


Pre-pressurization of the system increases the incipient 
superheat by dislodging or deactivating potential 
nucleation sites (flooding the cavities). 


Subcooling the refrigerant pool for a short period (<1 
hr) increases the incipient wall superheat to a value 
close to or even higher than that obtained with overnight 
settling time. Again this is due to deactivation of 
potential nucleation sites. 


Alkylbenzene oil tends to clog the re-entrant cavities of 
the HIGH FLUX tube, resulting in a higher incipient wall 
superheat. Cleaning the porous tubes with a brush to 
unclog the cavities returned the surface to its original 
state. 








VII. RECOMMENDATIONS 


Further studies should be conducted to verify the 
incipient behavior of 3% oil on the GEWA-K 19 and 26 fpi 
tubes and HIGH FLUX tube to try to see why incipient wall 
superheat was reduced. 


Conduct similar incipient tests on other remaining 
enhanced tubes to better understand the effect of oil 
addition on enhanced surfaces. 


Modify the existing system for multiple (2 and 3) tube 
bundle experiments and conduct tests to study the effects 
of tube arrangement position and orientation on 
incipience from other tubes. 


Conduct further tests on effect of subcooling, dissolved 
gases and pressurization on incipient wall superheat. 


Conduct studies of oil/R-124 mixture properties. 
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APPENDIX A: LIST OF DATA FILES 





epee | 
od 
[pata 
po203csMo eT 


[D0204sM0_ 
aby 








































(eee eer aeeenrere TEER 


















jrruz mame | ruse Tyre 
|oozi0asm3 | smooth 
| 0215sé3 
. 
pata 

Data | 
D0223BGK193 
D0223CK193 Data 















[ere snares 
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Data | 


| rxue name 
| p0224cK193 


DEGASSING 





D0224AGK193 
D0225GK193 


D0225BGK1910 
D0226GK1910 GEWA-K 19 


GEWA-K_19 
GEWA-K_19 
GEWA-K_19 


GEWA-K 26 
DO306AGK263__| GEWA-K 26 
GEWA-K 26 


poso7cK263 | czwa-x 26a | Data 


ALS 


Te OO, ee eee ee 


Hd 



































PURPOSE 
D0308AGK26: 
D0308GK2610 

Data 


Data 
D0311GK2610 























ee 


D0311BGK2610 
[Data 
D0313BHFO 


Dat 
Data 
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DEGASSING 


fpurrose | 
[bata 


Data 
Data 


FILE NAME 








[possismo_— | smoot [wf ata 
|possiasmo | smooms =f ata 
D0331BSM0 Data 


D0331CSM0 SMOOTH 


D0331DSM0 SMOOTH 
DO0401SM0 SMOOTH Data 


D0401asMo SMOOTH [Bi Data 





Data 
Data 


D0401BSM0 Data 
po401CSMO 
p0401DsMO 
p0404SMO 


SMOOTH Data 


Data 
Data 





Data 


t 


D0404ASM0 





a2 














FILE NAME TUBE TYPE DEGASSING 
D0405SMO0 Data | 
Do405AsMo 
D0405BSM0 
D0405CSMO SMOOTH Data 

D0405DSM0 SMOOTH lpata 
D0405ESMO SMOOTH 
D0405ASM3 SMOOTH 
D0406SM3 SMOOTH 
D0406ASM3 SMOOTH 
D0406BSM3 SMOOTH 
D0406CSM3 SMOOTH 
D0406DSM3 SMOOTH 
D0406ESM3 SMOOTH 
D0406FSM3 SMOOTH 
D0407SM3 SMOOTH 
D0407ASM3 SMOOTH 
D0407BSM3 SMOOTH 
D0407CSM3 SMOOTH 
D0407DSM3 SMOOTH 
D0408SM3 SMOOTH 
|ooaogasms | smoomm tw ata 
|ooacessms | smoomm_ fw ata 
[ooaoscsma__——| smoorm_ fw ata 
[soacenaa [amor —_[p |oste __| 
[poaogsus | smoorn_ fw ata 
fpoaogsmio | smoomm = fp fata 

















L22 




















































DEGASSING 


FILE NAME 





n 
{ 
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[rruz name | ruse Tyre _| DEGASSING 
GEWA-K 19 

GEWA-K_19 

| boa2opcKi910 _|cewa-k 19 | 8 
[boazonro | wicH-rLux |B 
D0421HFO 
D0421AHFO 
D0421CHFO 
D0421DHFO 
D0421EHFO 
D0422HFO 
D0422AHFO 


D0422BHFO HIGH-FLUX 


D0422CHFO 









Data 
pata 


Uncertainty 
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Ces Cae 
PBS ee ees 
Pe Si ata 
fp fata 
|poa23cur3— | nrcu-ruux |e | ata 
ee 
ce 
pe 4 
Cie! 








[ooazacurs | wrcu-ruux |e pata 
jooa2suri0 | wrcu-ruux [pata 
jooa2spuri0 | urcu-ruux |p aca 
jooazsanrio | urcu-ruux |e | pata 
pooascurio | urcn-roux [Bf bata 
jooazeurio so urcu-ruux |e ata 
jpoa2cauri0 | urcu-ruux fe | Data 
D0426BHF10 Data 
D0427GK260 
D0427GK260 
D0427AGK260 . 
D0428GK260 
D0428AGK260 
D0428BGK260 
D0428GK263 
D0428AGK2 63 
D0429GK263 





Data 
Data 


ll 
iii 


Data 
Data 
Data 
Data 





Data 






Data 


Data 





| 
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D0430GK2610 
D043 0AGK2610 


D0504GK190 





The tabulated data was filed using the same filing system used 
by Bertsch [Ref. 6]. For example filename ‘D0312AGK2610': 
1. First letter D refers to data. 
2. The next four numbers refer to the month and day. 


3. The next letter represents the order data was taken for 
a particular day; 


no letter - first run of the day 
A - Second run 
B - Third run 
etc 


4. The next two to four letters indicate tube type; 


SM = Smooth 

GK19 = GEWA-K 19 fpi 
GK26 = GEWA-K 26 fpi 
HF = HIGH FLUX 
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5. The last one to two numbers indicate the percentage of 


oil used; 
0 = 0% oil 
3 = 3% oil 
10 = 10% oil 


Thus, the example filename was a test data conducted on 
Mar 12, and was the second run for the day using a GEWA-K 26 


fpi tube with 10% oil. 
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APPENDIX B. REPRESENTATIVE DATA SET 


Oete : 30 Mer 1993 


NOTE: Program neme : ORPGB 
Oisk number = 0@ 
New file name: 0@33GESMe 
TC ars defective at lecations 7 a 
Tube Number: 4 


Oata Set Number = { @Buik O11 2 = 9.2 
TIME: 14:46:S$3 
TC No: t 2 3 4 s 6 7 8 
Teno : 4.38 4.46 4.45 4.481 4.5) 4.16 -99.99 -99.99 
Twa Tliad Thiad2 f[fvaer Psat Tsuno 
4.36 3.14 2.5 2.10 10.76 <-it.1 
Rho =) (1425.496 
Mu = 008352 
K = 07485 
Co = 1061.21519 
Tfhilm © 3.22636 


Hfg = 19§89S7.68568 
Tsat © 2.99544 
Hber 306 «160.71038 
Thetab Htube Qdo 
2.262 2.449€+92 S.SITE+O2 


Oata Set Number = 2 Bulk 011 % = 6.2 
TIME: V12tS:43 
TC No: I 2 3 4 Ss 6 7 8 
Temp : 5.00 S.t! §.14 S.@7 6.94 4,76 -99.99 -99.99 
Twa Tliad Tliqd2 Tvapr Psat Tsump 
4.99 2.18 2.2! 2.16 (@.82 <-tl.t 


Rho =) «(1 424.482 
Mu = .@003312 

K = .07382 

Ca = 1061.9836S 
Tftin * 3.97395 


Hfg # 1§8929.45867 
Tsat =~ 2.15584 
Hber = 1780.97572 
Thetab Htube Qdo 
2.836 3.2@SE+@2 9.091E+02 


Deta Set Number = 3 Bulk O11 2 = @.9 

TIME: 192t§:eg 

TC Ne: 1 2 3 4 S 6 7 8 
Temp : 7.13 7.34 7.38 =7.33 7.2@ 6.87 -99.99 -99.99 
Twa Tliad Tliqd2 Tvapr Psat Tsump 


7.19 2.22 2.25 2.17 10.83 <-11.1 
Rhea =) 1421.247 
Mu = .0003258 
K 2 07353 
Co © 1064,44863 
Tfilm = 4.68168 
Hfg = =_«-«158923.45200 


T3et 27 T8869 
Heer = 199.93379 
Thetab Htube Qdp 
S.0@26 3.24@E+@2 | .&629E+03 


128 


Date Set Number = 4 Bulk 011 % « @.@ 
TIME: 11:17:86 
TC No: 1 2 3 4 5 6 7 8 
Temp : 10.48 18.61 10.85 10.57 1@.22 9.66 -99.99 -99.99 
Tue Titad Tliad2 Tvapr Psat Tsump 

10.48 2.21 2.25 2.13 10.88 -11.2 


Rho = «=1416.588 
fu = .8083207 

K = .07296 

Cp = 1068.01306 
Tfilm = 6.26510 


Hfg = 1§8940.26987 

Tset © 2.1327) 

Hber = 229,63708 

Theteb Htube Qdp 

8.265 3.3S8E+@2 2.775E+93 


Reta Set Number = S Bult O11 2% = @.8 
TIME: 14:19:39 
TC No: 1 2 3 4 s 6 7 8 
Temp : 16.06 16.76 16.8% 16.76 16.33 15.41°-95.99 -99.99 
Twa Tliad Thiad2 Tvepr Pset Tsurp 
16.33 2.28 2.29 2.13 18.79 -11.2 


Rho == 1407.746 
Yu = .0803096 

K = .07191 

Co = 1074.81715 
Tfilm = 9.23281 


Hfg = 158940.87046 

Tsat = 2.13142 

Hbar = 267.36095 

Thetab Htube Qdp 
14,203 3.@77E+@2 4.371E+03 


Dete Set Number = 6 Bult O11 2 = 0.0 
TIME: 11:22:02 
TC No: 1 2 3 4 Ss 6 7 8 
Temp + 20.65 21.47 21.53 21.23 20.32 18.63 -99.99 -99,99 
Twe Tliad Tliad2 Tvapr Psat Tsump 
20.63 2.15 2.19 1.97 10.65 -11.2 


Rho = 14801.495 

Mu = .0003023 

K = 07117 

Cp = 1078.64687 
Tfilm = 11.30038 

Hfg = 1§9014.11805 
Tsat = 1.97460 

Hber © 283.07728 
Theteb Htube Qdp 


18.652 3.258&+02 6.077E+@3 
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Geta Set Number = 7 Quik 011 2 @ 8.2 


TIME: 11323:39 
TC No: ! 2 3 4 5 5 ? 8 
Temp : 25.2! 26.41 26.59 26.24 25.37 23.72 -99.99 -99.99 
Twe (liad Tliad2 Tveor Psat Tsump 

25.53 2.32 2.38 2.17 10.83 -11.3 

Rho ==) 1393.88@ 

Mu = ,0062935 

K = .07825 

Cp = 1085.69564 

Tfilm © 13.8496! 

Hfg = 1§8922.85131 

Tset = 2.16998 

Hoar = 306.62237 

Thetad Htube Qdo 

23.359 3.272E+@2 7.642E+03 
Data Set Number = 8 Bulk O11 x © @.@ 

TIME: 11:26:@S 

TC No: t 2 3 4 S =) 7 8 
Temp : 10.95 11.28 tt.er 10.68 1@.41 1@.53 -99.99 -99.99 
Twa Tliad Tliad2 Tvape Psat Tsure 

10.74 2.22 2.25 2.23 10.89 -11.3 

Rho == 1415.933 

Mu = .@@@3198 

« = 87283 

Ce = 1068.S155¢ 

Tftlm = §.48669 
Hfg = 1§8894.914658 

Tsat © 2.23166 

Heer = 231.54919 

Theteb Htube Qdo 


Oata Set Number = 16 Quik O11 2 = 0.9 


TIME: Ih: 3@24§ 
TC No: ! 2 3 4 Ss 6 7 8 
Temo : 13.12 13.60 13.34 12.61 12.26 12.44 -99.99 -99.99 
Twe Tliad Tliqad2 Tvaor Psat Tsureo 

12.78 2.28 2.33 2.28 10.93 -11.3 

Rho ==) 1412.84 

fu = .0003!1S9 

K = .07251 

Co ~ 18070.89183 

Thila = 7.52908 

Hfg © 1S6672.382ee 

Tset = 2.27791 

Hber = 245.59642 

Thetab Htube Qdp 

1@.5@2 1.807E+@3 1. a97E+04 
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Date Set Number = 11 Bulk O11 % = 0.2 
TIME: 11:36:29 
TC No: t 2 3 4 5 6 7 8 
Temp : 15.67 16.2@ 15.79 15.@2 14.69 15.11 -99.99 -99.99 
Twa Tliad Tliqd2 Tvapr Psat Tsump 
19.2) 2.2) 2.24 2.19 180.85 -1t1.2 


Rho ==) 1409.33! 
Mu = 0003116 

K = 07209 

Co = 1073.S9415 
TFilm = 8.70433 


Hfg ss =_-:158911.43766 

Tsat = 2.19439 

Hbar = 260.88268 

Thetab Htube Qdp 
13.02@ 2.S97E+03 3.382E+04 


Oata Set Number = {2 Bulk 011 2% 0.0 
TIME: 11:38:21 
TC No: 1 2 3 ‘4 5 6 7 8 
Temp : 18.04 18.31 18.06 17.47 17.19 17.85 -99.99 -99.99 
Twa Tlind Tliqd2 fTvapr Psat Tsump 
17.49 2.25 2.3) 2.26 «61@.9t -tt.? 


Rho = 1405.826 
Mu = .0803073 

K = .07168 

Cp =< 1076.29935 
Tfilm = 9.87057 


Hfg = = -1$8882.59790 

Tsat = 2.25607 

Hbar = 272.90109 

Thetab Htube Qdp 
15.229 3.828E€+63 5.829E+04 


Data Set Number = 13 Bulk O11 2 = 8.0 
TIME: tisats@4 
TC No: 1 2 3 4 5 6 7 8 
Temp : 19.57 19.56 19.48 16.96 18.74 19.51 -99.99 -99.99 
Twa Tliad Tliqd2 Tvapr Psat Tsump 
16.81 2.22 2.28 2.19 18.85 -18.9 


Rho = 1403.932 

fu = .8@0305! 

K = (8794S 

Cp = 1877.762S7 
Tfhilm = 10.4973! 

Hfg e 1§8915.0421@ 
Tsat = 2.18668 

Hbar = 279.769S2 
Thetab Htube Qdp 


16.621 S.1S6E+03 8.570E+04 
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APPENDIX C. SAMPLE CALCULATIONS 


Data set number 10 of Data file D0330ESMO (smooth tube 
with 0% oil concentration) was used for the sample 
calculation. 

A. TEST-TUBE DIMENSIONS 

D, = 0.015875 m 

D, = 0.0127 m 

D, = 0.012446 m 

Lie = 052032. 'm 

L, = 0.0762 m 
B. MEASURED PARAMETERS 


Vv 


V, x 25 = 4.444 x 25 = 203.25 volts 
I 


I, x 1.9182 = 0.9536 x 1.9182 = 1.829 amps 
T1 = 23.12 °C 

T2 = 13.60 °C 

T3 = 13.34 °C 

T4 = 12.61 °C 

T5 = 12.24 

T6 = 12.44 °C 

Tsat = 2.28 °C 


kc = 344 W/mK 
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C. OUTER WALL TEMPERATURE OF THE BOILING TUBE 
1. Surface Area of Smooth Tube (Ac): 
p=" x Do = ® x .015875 m =.04987 m 


Ac = ®(Do*-Di*)/4 = m((.015875)?-(.0127)-) /4 


Ac 7.13x10°° m 
2. Heater Power (Q,): 
Q, = VI = 111.15 volts x 1.829 amps = 203.25 W 
3. Tube inside Temperature (Tavg): 
Tava = Yr sn = Yr/6 = 12.89 °C 
4. Tube outside wall temperature (T,,): 

T,, 18 calculated by knowing the inside wall 
temperature using Fourier’s Conduction Law, assuming uniform 
radial conduction. 

Ty = Tug = [ oh: (1n(Do/D1) / (2m L- ke) ) 


Tig = 12.89: = [ 203.25- (1n(.015875/.012446) / (2m L- kc) )] 


Tey = 12°89 = 2.11 = 12778 °C 
5. Wall Superheat (aT): 


aT =(T,.- Tyr) = 12.78 - 2.28 = 10.5 


D. PROPERTIES OF R-124 AT FILM TEMPERATURE: 

All thermophysical properties of R-124 were computed at 
film temperature using the REFPROP equations as provided by 
Bertsch [Ref. 6]. 

1. Film Temperature (T,): 


Tf£= (Two+Tsat) /2=(12.78+2.28)/2= 7.53 °C 
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Dynamic Viscosity (1): 
w= 3.4593~-(.0426) T,+(3.9485x10-4) T,?- (4.193x10°°) T,?+ 
(2.0709x10-°) T,4 
Substituting T;: 
pt = 3.4593-(.0426)-7.53+(3.9485x10)- (7.53) °~ 
(4.193x10°°)- (7.53)?+ (2.0709x10°*)- (7.53) * 
wt = 3.159x10* N-s/m’ 


Density (p): 


p = 1434.8-(2.8619)T,-(6.7267x10-?) T,?- (7.2852x10°) T,’ 
p = 1434.8-(2.8619)-7.53 - (6.7267x10°°): (7.53)? 

- (7.2852x107°)- (9.2)? 
p = 1412.837 kg/m’ 


Kinematic Viscosity (v): 

D = H/p = 3.159x104/1412.837 = 2.236x10°’ m/s 

Thermal Conductivity (k): 

k = 7.5191x107? - (3.5436x1074).T, - (1.9545x10°7) T;,? 
* (3.1835x10"") 2TY 


k = 7.5191x10-?- (3.5436x1074)- (7.53) 
- (1.9545x10°7)- (7.53)?+ (3.1835x10°°)- (7.53)? 
k = 7.251x10°? W/m: K 


Specific Heat (C,): 

C, = 1.0542+(2.1405x10™°) T£+(1.0709x10~°) TE? 
~(6.4721x107°) TE- (1.4324x10°°) TE4+ (4.136x10°'') TE? 

Cp = 1.0542+(2.1405x107*)- (7.53) 
#° (1.070910) "(7-532 )*=(6.4721%10"")* (922) * 
-(1.4324x107°)- (9.2) 4+(4.136x107")- (9.2)5 
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C, = 1070.89 J/kg: K 
7. #‘YThermal Diffusivity (a): 


a 


k/(p-C,)= 7.251x10*/ (1412.84 -1070.89) 


a 4.792x10°° m?/s 


8. Volumetric Thermal Expansion Coefficient (f): 


B = ~ (Ap/AT) /p = (Puc=Peac) 7 (Poc=Taae) 
= -(1396.97 - 1428.24)/(12.78 - 2.28) 
B = 2.107x10° (1/K) 


9. Prandtl Number (Pr): 


Pr = vV/Q@ = 4.666 


E. HEAT FLUX (Q,) CALCULATION: 

To find the average natural convection heat-transfer 
coefficient at the non-boiling ends of the test tube, the 
Churchill and Chu [Ref. 45] correlation was used: 

For simplification let R represent the following in the 
correlation : 

R=[ (g: B: Do?-©@- tanh (m: Lu) ) / (8 a: Lu: m) }!/° 
where: m = [(h'p)/(kce-Ac) ]? 
then: 

h = k/Dol .6+.387({R/(1+(0.559/Pr) 6) 9/27}2)| 
The value of h was computed iteratively. An initial value for 
h of 190 (W/m’-K) was used. The following results were 
obtained: 


h 245.59 (W/m? K) 


m = 24.93 (1/m) 
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therefore: 








Q, = (h:-p:kce-Ac)!/?-e-tanh(m.Lu) = 5.50 W 
HEAT FLUX THROUGH ACTIVE BOILING LENGTH 
Q=QQ, - 2° = 203.25 - (2°5.50) = 192.25 W 
A, = ™-D,-L = m-0.015875-0.2032=1.0134x107 m 
q" = Q/A, = 192.25/1.0134x10° = 18970.79 W/m? 
h = q"/O = 18970.79/10.5 = 1806.74 W/m*-K 


The following results were obtained from the computer 


using the data acquisition and reduction program DRPGP: 





q" = 18970 W/m’ 
e 
h 


102502) °C 


1807 W/m?*- °K 
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APPENDIX D. UNCERTAINTY ANALYSES 


The same data run used for Appendix C for a heat flux of 
18970 w/m’ was used for uncertainty analyses, thus all the 
measured and most of the calculated results still apply. A 
higher heat flux of 85700 w/m’ (data set number 13) was also 
analyzed to show the difference between the two heat fluxes. 
All uncertainties are presented as a percentage of the 
calculated parameter, i.e. a relative uncertainty. The 


uncertainty analysis method suggested by Kline and McClintock 











{(Ref. 47] was used. for example: R = R(xX,,X%2,X%3,-- ».-X,), 
then: 
8R,,\ ,(8Ra,\ rx, \']2 
8R = (gee%,] “(¥ 03, wee + (gEox| 
where: 


dR = uncertainty of the desired dependent variable 
X, = measured variable 


6x, = uncertainty in the measured variable 


A. HEAT TRANSFER RATE UNCERTAINTY: 


0.025 amps 


i; 0.9536 amps bi, 


4.444 volts dv, 0.05 volts 


S 
" 
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so we have: 

















50, - [ 0.05 |? 


B. SURFACE AREA UNCERTAINTY: 
The following dimensions were 
manufacturer’s data sheet, 


D, = 15.88 











(mm) 5D,= 0.1 (mm) 
L = 203.2 (mm) dL = 0.1 (mm) 
A, = Do x L 
1 
BA, _ |(5D.)° , /8L\?|? 
A, D, Ga 
8A, 0.1 \2 0.1 \2]3 
A, (Gea + (so375) | ieee 
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east acces oi SAR lea ee St need i indie iM sac ni te SEE Stan anda a dct SSE Sac lS en Spa ire ESS 


1 
0.025 \2]3 
4.444) * Sees) | a eana® 


taken 





from the 


and from actual measurements: 





Cc. 


Note: 


WALL SUPERHEAT UNCERTAINTY 


The following temperatures were measured: 


T, = 13.12 °C 
T, = 13.60 °C 
T; = 13.34 °C 
T, = 12.61 © 
tT; = 12.24 °c 
T, = 12.44 °c 


Tage 12.89 °C 
Standard Deviation = 0.541 °C 


Ter = 2.28 °C and 8T,,, = 0.01 °C 
D 
1n(=*) 
T= Tag - 24 \—— 
wo Sa MRO Kk, 


The second term on the right hand side of the above 


equation is the Fourier heat-transfer conduction term, and 


can be neglected for the uncertainty analysis since its value 


is much less compared to the temperature standard deviation, 


SO; 


Te oT S02 69° °C 
AT = T,, - Ty = 12.89 °C - 2.28 = 10.61 © 
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tole 

















substituting values: 


BAT _ [/.0.541\? (3 01 : 
AT (Grorst 61 (a0 ool eae 


D. HEAT FLUX UNCERTAINTY: 


(QO, — 2 X O-) 


q* A, 


where: Q, =Vx I_ = 203.25 watts 


and 50, = (203.25) (0.0285) = 5.793 watts 
q = 18970 w/m (measured) 
A, = 0.010137 m 


Qn ~ 2.Q¢ = q.A, = (18970) (0.010137) = 192.29 W 
solving for Q;: 

Q, = (203.25 - 192.29)/2 = 5.48 W 
assuming the same proportion uncertainty for Q,, we have: 


50, = (5.48) (0.0285) = 0.1562 W 


and now for the uncertainty, we have: 


7 80, \ (2-80 \ | (8a) 
O, ~— 2-0¢ On ~ 2-0¢ A, 


$¢ 
q 

















and substituting the computed values: 


1 
2 —_ 
2 3 
(372 2, (2.x 0.1562\? + (6:38 x20 *} | = 3.0% 


8g . 
gd 192.29 192.29 0.010137 





E. HEAT-TRANSFER COEFFICIENT UNCERTAINTY: 
h = q/DT 
and 


eo (S) aelT 


and substituting computed values we have: 


3h 3 


a (0.003)? + (0.005)7]2_ = 5.11% 


The above series of calculation was also conducted on the 
higher heat flux for the same data filename, furthermore an 
uncertainty analyses was also conducted on a high flux run 
with 10% oil at both low and high flux, the results are 


tabulated in Table C.1. 


Table C.1 RESULTS OF UNCERTAINTY ANALYSES 


| 
| 


i 
1 
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